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SUMMARY 


A facility has been developed for dynamic testing of straight and 
convergent-tapered seals, with the capability of separately determining both 
direct and cross-coupled stiffness, damping, and added mass coefficients. The 
test apparatus causes the seal journal to execute small-eccentricity centered 
circular orbits within its bearings. Dynamic measurements are made and 
recorded of the seal-displacement-vector components, and of the pressure field. 
The pressure field is integrated to yield seal reaction-force components. The' 
displacement and force vector components are analyzed via a generalized 
Newton-Raphson procedure to yield the desired seal dynamic coefficients. 
Representative test data are provided and discussed. 


INTRODUCTION 


Black [1, 2, 3, 4] in a series of publications incorporating theoretical 
and experimental results has demonstrated that the rotordynamlc behavior of 
pumps IS critically dependent on forces developed by neck-rings and Interstage 

in figure 1. Subsequent experience [5] has demonstrated 
that the stability of cryogenic turbopumps is comparably dependent on seal 
orces. The test program discussed here was stimulated by stability diffi- 
culties encountered in developing the turbopump of [5], and has the objective 

of separately identifying stiffness, damping, and Inertia coefficients for 
turbulent seals. 

The contents of this section are provided to briefly review theoretical 
models and prior experimental results and procedures for seals. Also, since 
t e seals of figure 1 are geometrically similar to plain journal bearings, 
applicable prior test programs to identify journal bearing coefficients are 


*The work reported herein was supported by NASA Lewis under NASA Grant 3200- 
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Seal Analysis: Leakage and Dynamic Coefficients 

Black [1, 2, 3, 4] is largely responsible for developing currently 
employed dynamic seal models. Black's analysis yields a definition of the 
force acting on a rotor due to its motion at a seal location, and is based on 
the following leakage relationship from Yamada [6] for flow between concentric 
rotating cylinders. 


AP = (1 + ? + 2a)pV2/2 


( 1 ) 


where C is a constant entry- loss coefficient, p is the fluid density, V is the 
average fluid velocity, and o is a friction-loss coefficient defined by 


a = AL/C 


( 2 ) 


In the above, L is the seal length, C is the radial clearance, and A has been 
defined by Yamada to be the following function of the axial and radial Reynolds 

numbers (R , R ) 
a r 


A = 0.079 R 


a 


; R = 2VC/v, R = Ro)C/v 
a r 


(3) 


where v is the fluid's kinematic viscosity, R is the seal radius, and o) is the 
rotor's rotational speed. The friction law definition of Eq. (3), Yamada s 
definition for X, is based on an assumed 1/7 power velocity distribution, and 
fits the Blasius equation for pipe friction. 

Black's analysis for a plain non-serrated seal yields a motion/reaction- 
force definition of the form 


A 

RttAP 


R 



K k r C c r m o r 

[,-']{)+[- 7-H. } + [- -H- } 

-k K r^ -c C r^ o m r^ 
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where 


K = 




2x2/4, ic = 


y ^ oaT / 2 , 


C = 




c = 




m = 




T = L/V (5) 


In Black’s original analysis [1], the coefficients p^, ho 

developed for short seals for which, "circumferential pressure-induced flows 
are negligible compared with axial flows". Black subsequently [2] examined 
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the effect of circumferential pressure-induced flow for finite-length seals, 
and developed the following formulae to account for finite (L/R) ratios 


+ 0.28(L/R)2}-i'' 
Pl(|) = + 0.23(L/R)2}-1 > 

P2^-|) = ^2^^ 0.06(L/R)2}-1^ 


Bla^k's^sec^nd refinement of the original development [3] was the defini- 
tion of uqj Pj, P2 iri terms of the following additional parameter 


3 


7R 7R 

a a 


(7) 


which accounts for a circumferential variation in A due to a radial displace- 
ment perturbation from a centered position. Plots of pQ, pj, and ^2 ^re pro- 
vided in figure 2 as a function of 3 and a for C = 0.5, These coefficients 
are comparatively insensitive to anticipated variations of the entrance loss 
factor 5. 

Finally, Black [4] examined the influence of inlet swirl on seal coeffi- 
cients. Specifically, in previous analyses, a fluid element entering a seal 
was assumed to instantaneously achieve the half-speed tangential velocity 
Rm/2. Black in [4] demonstrates that a fluid element must travel a substantial 
distance axially along the seal before asymptotically approaching this limiting 
velocity. The practical consequence of this swirl effect is that predictions 
for the cross-coupling terms k, c may be substantially reduced. 


One of the authors [7] has recently completed a seal analysis based on 
Hirs turbulent lubrication model [8,9], which largely repeats Black’s develop- 
ments (which were based on various ad hoc models). The results resemble, but 
do not coincide with Black's results. They also do not include the finite- 
length correction of Eq. (6). 

Prior Seal Testing Procedures and Results 


The pertinent data which must be measured to confirm the seal leakage 
model of Eqs. (1) - (3) are AP , V (from flow rate), m , and the axial pressure 
gradient within the seal. This latter measurement yields a which in turn 
yields X. Yamada's model for the friction factor was based on testing for 
these variables over the Reynolds number range (200 Ra £ 40,000; 0 < R^- < 
40,000) and clearance to radius ratios of (.0106 £ C/R < Toi 29 ) . ~ ~ 

Various approaches can be taken to the measurement of seal dynamic pro- 
perties as defined by Eqs. (4) and (5), For example, if the journal segment 
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of the seal is stationary (i.e., - fy - - ry 0), Eq. (4) can be 

inverted to obtain 


{ 



X 

RttPK 2 
e 


K 



R 



K 2 
e 




Hence by applying the static load definition (F^ = , Fy = 0), and measuring 

the displacement components r^, Ty, one obtains a combined measure of the 
direct and cross-coupled stiffness coefficients. This is predominantly the 
type of testing performed by Black, who cites results in the form of 
"receptance magnitudes", i.e.. 


Ir |/F 
' ' s 


1/2 

X{K^ + k^l /RttAP 


From Black's model, the relative magnitudes of the direct K and cross-coupled 
R stiffness coefficients depend on the relative magnitudes of the axial and 
radial Reynolds numbers. Specifically, at zero running speeds k is zero, but 
increases with R;j- , and can exceed K. 

For softly supported rotors, the direct stiffness of a non-serrated seal 
represented by K may significantly influence the location of a critical speed. 
However, the stability of a flexible rotor is less sensitive to the direct 
stiffness term, depending primarily on the cross— coupled stiffness coefficient 
k, and direct damping coefficient C. For most rotors, the cros.s— coupled 
damping coefficients c, and inertia terms tfi have no appreciable influence on 
rotor stability or response. 

Most of Black’s testing [2, 3, 10, 11] has been of the static nature cited 
above. The second type of test cited consists of analytically modeling a test 
rotor including the theoretically predicted seal dynamics, and comparing the 
dynamic characteristics of the model with test data. For example, in [2] the 
test rotor was rapped and a correlation was made with the observed logarithmic 
decrement on the decay curve. In [3], known imbalances were applied to the 
test rotor, and a comparison was made with synchronous amplitudes and phase, 
cr i t ical” speed location, and onset speed of instability. Comparisons between 
rotor model results and tests, of this nature, are helpful in deciding whether 
the general seal model is reasonable. However, this type of test— correlation 
does not yield specific information about the individual dynamic coefficients. 
Further, discrepancies between predictions and results can be the result of 
either an inadequate rotor model or an inadequate seal model. For example, 
in [3] , Black indicates that discrepancies in synchronous amplitude and phase 
results could result from an inadequate initial balance. 

A summary of the test results of references [2, 3, 10, 11] is provided in 
Table 1. The correlation in these tests ranges from "good" to "fair". The 
nature and results of the test support the following general conclusions 
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concerning the adequacy of Black's dynamic seal model: 

(a) Over the Reynolds number range tested, the prediction of the direct 
stiffness coefficient K is adequate for plain and serrated seals, 
although less accurate for serrated seals. Black's test results 
indicate a divergence between tests and theory for the direct damp- 
ing coefficient C as the axial Reynolds number is increased. 

(b) Although the data cited generally supports Black's dynamic seal 
model over the Reynolds number range considered, it is inadequate 
to specifically verify the proposed relationships [Eq. (5)] for the 
dynamic coefficients as functions of the axial and radial Reynolds 
numbers . 

Identification of the dynamic coefficients of seals in a centered position as 
functions of the axial and radial Reynolds numbers is the objective of the 
current test program . 

Prior Journal-Bearing Coefficient Identification Approaches 

The seals of figure 1 are geometrically similar to plain journal bearings 
but have larger C/R ratios on the order of 0.01 as compared to bearings for 
which C/R is on the order of 0.001. Seals customarily operate in the turbulent 
regime, both axially and circumferentially, and have a substantial direct 
stiffness at a centered zero eccentricity position. Further, seals are 
nominally designed to operate in a centered position, while the operating 
eccentricities of journal bearings vary with running speed and load. Hence, 
dynamic bearing-identification work has generally had the objective of validat- 
ing dynamic coefficients versus eccentricity relationships. The general 
similarities between bearings and seals are such that procedures developed 
for bearing-coefficient-identification may also apply for seals and are 
briefly reviewed below. 

On the basis of various analyses, the motion/reaction-force relationship 
for a hydrodynamic bearing is defined, for small motion about an equilibrium 
position, by the equation 


C C 
XX xy 
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( 8 ) 


The equations 
two identical 


of motion of a rigid rotor of mass 2M supported symmetrically by 
bearings can then be stated 
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Where F^, Fy are the components of the external force vector, a is the 
imbalance vector magnitude and w is the constant rotor spin speed. 

Published work related to the identification of the stiffness and damping 
coefficients of Eq. (9) date from Gleinecke [12] who excited the bearing seg- 
ment of a 120 mm model bearing in two mutually orthogonal directions while 
measuring the amplitude and phase of the relative motion between the bearing 
and journal. The stiffness and damping coefficients are calculated from the 
frequency-domain equations. Morton [13] adopted this test procedure on a 
full-scale 308 mm (20 in.) industrial bearing, to calculate bearing stiffness 
and damping coefficients, and subsequently [14] developed a technique for 
introducing a step input into a full-scale operating turbine bearing, thereby 
defining the bearing transfer function. Burrows and Stanway [15] have proposed 
the use of a pseudo-random-binary sequency (pros) excitation force with a 
multiple-regression analysis for estimating the coefficients. Their estimation 
procedure consists of the following steps: 

(a) The governing differential equations of motion (2) are expressed 
in state-variable format, i.e., as a system of first-order 
differential equations. 

(b) The first-order differential equations are replaced by first-order 
difference equations. 

(c) The unknown coefficients are calculated based on a minimum error- 
squared criterion from measurements of the state variables and the 
pros input signal. 

They also applied the method to the estimation of squeeze-film damper coeffi- 
cients [16]. 


THE UNIVERSITY OF LOUISVILLE DYNAMIC SEAL TEST PROGRAM 


Test Section Design 

Figure 3 illustrates the test-section design employed in the current seal 
test program. Water enters the center of the section and flows axially across 
the two rotating test seals exiting at the bottom of the test section. The 
seal journals (L = 4 in = 10.16 cm, D = 2 in = 5.08 cm) are mounted eccentri- 
cally on the shaft with a constant eccentricity eg = .005 in = 1.27 x 10“^m. 
The nominal seal clearance is C = .02 in = 5.08 mm, which yields C/R = .01. 
Accordingly, shaft rotation causes the seal journals to execute circular 
centered orbits at the nominal eccentricity ratio e = 0.25. Axial and radial 
Reynolds numbers may be specified over the range R3s[5000, 30,000] 

R^e[0, 11,000] by varying the shaft rotational speed (0 - 4,000 rpm) and flow- 
rate. Shaft-speed is measured by a once-per-revolution counter, while turbine 
flowmeters separately measure flowrate through each seal. 

The rotor of figure 3 is supported in Torrington hollow roller 
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bearings [17]. These precision bearings are preloaded radially, have zero 
internal clearances, and an accurately predictable radial stiffness. When 
supported in these bearings, the stiff rotor design of figure 3 yields a first 
critical speed of approximately 12,500 rpm. The end thrust bearing is pro- 
vided to react the small axial load developed by the opposed test-seal design. 

Instrumentation and Data Analysis 

The dynamic instrumentation illustrated in figures 3 and 4 consists, for 
each seal, of Bently eddy-current motion transducers and five piezo-electric 
pressure transducers which are distributed both axially and circumferentially 
a ong and around the seal. The circumferential "clocking" of the pressure 
yansducers is provided primarily as a matter of convenience, since the trans- 
ducers are provided exclusively to define the time history of the axial pres- 
sure distribution. Since the seal journal is forced to execute a closed 
circuyr orMt at constant speed m within its journal, the steady-state pres- 
suy ystribution is constant with respect to an observer fixed to the shaft 
and the circumferential pressure distribution at time T, p(s,6) is definable’ 

in yrms of either past or future time measurements p(z,t) at a fixed value 
Ox 0 , 


The direct extraction of circumferential pressure distributions from 
pressure time histories also permits the "reconstruction" of a pressure time 
ysyry at a g^en circumferential location (e.g. , 0 = 0), despite the fact 
yat the trayducers are distributed circumferentially around the seal as 
illustrayd in figure 4. Specifically, given the measurements pi(t), P 9 (t) 
••• P 5 Ct;, the corresponding pressure signals for 0 = 0 are ^ 


pf(t) = pgt) 

Pf(t) = p^Ct + ^/(O) 


pyt) = p^(t + 46/m) 


( 10 ) 


whey m IS ayy the shaft rotational speed. In words, the pressure measure 
ments p^(t) define the axial pressure distribution at time t = T for 0=0. 

The seal reaction forces at a given time t = T are defined by the 
integrals 


These bearings were donated by Tor ring ton through the kindness of W. L. 
Bowen, whose assistance is gratefully acknowledged. 
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_ 2tt L 

R (t) = - / / p(0,z) sin6 RdO dz 

^ 0 0 

_ 2tt L 

R (t) = - / / p(9,z) COS0 Rd0 dz 

y 0 0 


with 0 as illustrated in figure 4. These integrals may also be stated as 


2tt 


2tt 

0) 


R (t) = -RL / sin0 p(0)d = -RuL / sinwx p(x)dT 

^0 0 


y (12) 




2tt 

2tt _ _ 

-RL / COS0 p(0)d0 = -RwL / coswx p(x)dx 
0 0 


where p is the average axial pressure defined by 


p(0) = T ^ p(0,z)dz , p(t) = 7 / p(t,z)dz 


(13) 


The integral of Eq. (13) is evaluated numerically from the pressure time 
histories p*(t). The axial spacing of the pressure transducers has been chosen 
from Gauss-Legendre quadrature formulas [18] to minimize the error involved in 
evaluating this integral, and yields the following simple formula for the 
integration with respect to z 

p(t) = |{A^p*(t) + A^P*(t) + A^p*(t) + A^P*(t) + A^p*(t)} 


where 


= 0.23693, A^ = 0.47862, = 0.56889 


The time integrals in Eq . (12) are executed with a simple Simpson s-rule 

based algorithm. Note that a pressure time-history is required over the 
interval \t, t + 2tt/u] to obtain R^lTt) , RyOt)- Hence, pressure time histories 
p*(t) over a total time period AT will yield force histories Ry(t) , Ry(t) over 
the reduced interval [AT - 2 tt/w] . However, since the signals are periodic, an 
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adequate sampling rate^ yields sufficient output (rj^(t), ry(t), Ry(t) , RyCt)) 
from a limited number of data cycles to identify the dynamic seal coefficients 
K, k, C, c, m of the following dimensional verion of Eq . (4) 


m 0 
. 0 m 


C c 
“C C 




-k k J , r^ 


I R. 




X 


(14) 


Note Lnat m in this equation is the seal added fluid mass, and has nothing to 
do with the actual mass of the test rotor. Eq . (14) is solely the seal govern- 
ing equation, and the test procedure which yields rx(t) , ryCt), Rx(t) , RyCt) 
is not influenced by dynamics of the test-section rotor. 

Identification Procedure 

A generalized NR (Newton-Raphson) procedure [19] , [20] , [21] is employed 
for parameter identification. This procedure can be visualized as fitting the 
solution of Eqs. (14) to the observed data, in much the same manner regression 
analysis does with algebraic models. These solutions are obtained by numerical 
integration of the ordinary differential Eqs. (14), Note that this procedure 
operates on the differential Eq . (14), and is not a frequency-response or 
transfer-function approach. In fact, a frequency-response approach can provide 
only a restricted amount of information from the circular-orbit data of this 
program. This statement is illustrated by substituting the assumed solution 


r = A coscot , r = A sincot 

A Y 


into Eq , (14), solving for (k^, Ry) > and forming the dot and cross products 


*^xR| _ , ^ r-R 9 

- k - Cco , ~ - cco - K 


For a given speed (^, a frequency-response approach yields the sums on the right 
hand side of these expressions; however, the N-R procedure can separately 
identify the coefficients. 

The N— R procedure as applied to Eq. (14) requires the following first- 
order restatement 


Biomation data acquisition units are employed with a sampling rate per channel 
of 100,000 Hz to digitally record r^(t), ry(t) and Pj^(t). 
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Yl = Y3 




■ y/ 


^3 = 
^4 = 


-^5 >^3 - ^6 ^ ■ ^7 ^1 


+ye >'3 ■ ^'5 ^ ^ ^8 ^1 


^6 = ^7 = ^8 = ^9 = ° 


yg y2 + ^9 y 


h ^2 ^9 ^ 2 ^"^ 


(15) 


The identities relating Eqs . (14) and (15) are y^ = fy, Y2 = ^Y » Y3 " 

Y 4 = ry, Y5 = C/m, y^ = c/m, yy = K/m, yg = k/m, and y^ = 1/m. Note that the 
original differential equations have been augmented by the trivial differential 
equations, m = 0, K = 0, etc., to enforce the time invariance of these seal 
coefficients. The identification procedure ninimizes, in a least-square sense, 
^he eri“ors between the solution to Eq. (15) and the test data. 

Unfortunately, this procedure involves the solution of a multipoint- 
boundary value problem for which existence and uniqueness theorems are simply 
not available. Hence, the procedure was validated [21] by generating fake 
data with simulated noise and theoretically predicted seal coefficients, and 
then analyzing the data to "identify” the coefficients. Even with severely 
degraded data, the procedure yields errors less than 8%. 


TEST RESULTS AND DISCUSSION 


Experimental Results 

The beginning phase of testing has consisted of marginal advances 
followed by precipitous retreats. Tests were conducted on straight seals 
(D = 4 in, L/D = 0.5, C/R = 0.01) in June and October 1979. Because of data 
transfer problems, the bulk of the dynamic data taken in June was unusable, 
however, test results for three good R^, R^ sets are provided in the first 
three rows of Table 2. These data sets indicate that experimental estimates 
of m, C, c, k are smaller than predicted, while K is substantially larger. 
Figure 5(a) illustrates the theoretical [7] and experimental pressure distribu- 
tions for case 1 of Table 2, demonstrating that the experimental pressure 
amplitudes are smaller than predicted. The oscillation observable in Pg of 
this figure is exactly 18 times rotational speed (3660 rpm) , and is at present 
unexplained. However, as illustrated in figure 6(a), the integration which 
yields the force components eliminates this and all other fourier components, 
leaving only the fundamental component of running speed. The displacement 
vector components are illustrated in figure 6(b). The "notch" in the lower 
portion of these signals is the result of damage to the seal journal extension 
caused by rubbing of the displacement transducer probe during start up. 
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ho October 1979 were repeated in the straight seal configuration; 

however, the measured pressure fields were noisy and erratic, and one of the 
pressure transducers failed. When the test unit was disassembled, chemical 

bearing and Journal. These deposits apparently 
med when the test section was drained following the June test series. The 
deposits were irregular, and substantially Increased the surface roughness of 
e sea s. The results of two tests from this series are provided in the last 
two rows of Table 2, demonstrating the same basic trend as the June tests 
Note, however, that the added mass term is now larger than predicted. 


The leakage -AP data for a large number of tests 
deviated from Black's and Yamada's results in that the 
(by 5 to 8%) smaller than predicted. 


has consistently 
leakage is consistent 


ly 


Discussion of Results 

Generally speaking, an inadequate amount of data has been taken to 
support any conclusions about the adequacy of the theory. The differences 
etween the tests in June and October are consistent with the theory in that 
an increase in surface roughness would increase a. Over the range of a 
anticipated, the coefficient p2 is much more sensitive to a than pQ* PI 
Hence, a relatively sharper increase in the added mass m is anticipated 'with 
increasing surface roughness, 

results obtained in the tests as compared to Yamada's 
an lack s predictions are at present unexplained. It is possible that the 
oscillating entry clearance accounts for this result. 


Additional Planned Testing 

Tests were conducted on a seal with a convergent tapered sleeve segment 
during April 1980, and the dynamic data related to this test are currently 
being processed and analyzed. The taper angle for the seal tested is less 
t an optimal [22] from a direct stiffness viewpoint. Tests will also be con- 

.f """ optimal taper seal and a seal with a taper angle approxi- 

mately 25% larger than optimal. 


Tests will be repeated in June of 1980 on the straight seal over a wide 
Reynolds number range. The straight seal will then be modified by rounding 
the sleeve inlet. ® 


Strain-gauge pressure transducers are on order to replace and/or comple- 
ment the piezo-electric transducers presently employed. With these new trans- 
ducers, we should get a better idea of the dynamic entry loss, and dynamic 
pressure gradient. 
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TABLE 1. PRIOR DYNAMIC SEAL TESTS 


K k 


R - 18,317 

Theory 

47. 

14. 

9712 

3585 

0.06 

R^ - 10,079 

Estimated 

21. 

2.5 

14170 

2700 

0.013 

R - 11,997 

Theory 

31. 

14. 

4592 

3046 

0.06 

- 10,410 

Estimated 

22. 

0.33 

8859 

3204 

0.017 

R - 20,814 

Theory 

53.4 

14. 

12203 

3879 

0.06 

a 

- 10,079 

Estimated 

25.9 

2.89 

16972 

2556 

0.013 

R - 17,834 

Theory 

50.2 

4.0 

9106 

855 

0.07 

a 

- 2,644 

Estimated 

13.2 

14.3 

9662 

655 

0.098 

R - 19,022 

Theory 

54.0 

8.9 

10739 

2096 

0.07 

a 

R - 5.783 

Estimated 

18.6 

1.33 

14759 

2200 

0.113 


r 


Table 2. Test Results for June 1979 (Rows 1, 2. 3) and 
October 1979 (Columns 4, 5) versus Theory 17 J 
(ln-lb“Sec units) 
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